Relativistic stars with purely toroidal magnetic fields 
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We investigate the effects of the purely toroidal magnetic field on the equilibrium structures 
of the relativistic stars. The master equations for obtaining equilibrium solutions of relativistic 
rotating stars containing purely toroidal magnetic fields are derived for the first time. To solve 
these master equations numerically, we extend the Cook-Shapiro- Teukolsky scheme for calculating 
relativistic rotating stars containing no magnetic field to incorporate the effects of the purely toroidal 
magnetic fields. By using the numerical scheme, we then calculate a large number of the equilibrium 
configurations for a particular distribution of the magnetic field in order to explore the equilibrium 
properties. We also construct the equilibrium sequences of the constant baryon mass and/or the 
constant magnetic flux, which model the evolution of an isolated neutron star as it loses angular 
momentum via the gravitational waves. Important properties of the equilibrium configurations of 
the magnetized stars obtained in this study are summarized as follows ; (f) For the non-rotating 
stars, the matter distribution of the stars is prolately distorted due to the toroidal magnetic fields. 
(2) For the rapidly rotating stars, the shape of the stellar surface becomes oblate because of the 
centrifugal force. But, the matter distribution deep inside the star is sufficiently prolate for the 
mean matter distribution of the star to be prolate. (3) The stronger toroidal magnetic fields lead to 
the mass-shedding of the stars at the lower angular velocity. (4) For some equilibrium sequences of 
the constant baryon mass and magnetic flux, the stars can spin up as they lose angular momentum. 

PACS numbers: 



I. INTRODUCTION 



There is growing evidence for the existence of so-called magnetars, supermagnetized neutron stars with the magnetic 
10 14 - 10 15 G (e.g., 



fields of 
canonical neutron stars with 



23, 



42]). Although such stars are estimated to be only a subclass (~ 10%) of the 
10 12 — 10 13 G 22J, much attention has been paid to the objects because they pose many 



astrophysically exciting but unresolved problems. Giant flaring activities in the soft gamma repeaters observed for the 
last two decades have given us good opportunities to study the coupling of the interior to the magnetospheric structures 
of magnetars 3J, [35| , but the relationship between the crustal fraction and the subsequent star quarks has not yet 



been clarified (see references in HQ). The origin of the large magnetic field is also a big problem, whether generated 
at post-collapse in the rapidly rotating neutron star ^\ or descended from the main sequence stars 9]. Assuming 
large magnetic fields before core-collapse, researches have recently performed extensive magnetohydrodynamic (MHD) 



stellar collapse simulations 



0QHE 



31 



43j in order to understand the formation mechanism of magnetars, with 



different levels of sophistication in the treatment of equations of state (EOSs), the neutrino transport, and general 
relativity (see 2l[ for a review). Here it is worth mentioning that the gravitational waves from magnetars could be a 
useful source of information about magnetar interiors [^Q. From a microscopic point of view, the effects of magnetic 
field larger than the so-called QED limit of ~ 4.4 x 10 13 G, on the EOSs (e.g., 23j ) and the radiation processes have 
been also extensively investigated (see [3] for a review). For the understanding of the formation and evolution of the 
magnetars, it is indispensable to unify these macroscopic and microscopic studies, albeit not an easy job. 
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When one tries to study the above interesting issues, the construction of an equilibrium configuration of magnetars 
may be one of the most fundamental problems. In spite of the many magnetized star investigations that have so far 
been carried out, with different levels of sophistication in the treatment of the magnetic field structure, equations of 



state (EOSs), and the general relativity 



els oi sopmstication m tne treatn 
ivity 16, 17, IS, 36,0,Q 



36, 39, 40], there still remains room for more sophistication 



in the previous studies in which stellar structures are obtained within the framework of Newtonian dynamics. Since 
we are interested in magnetar structures, the general relativistic treatment due to the strong gravity of magnetars 
is required in constructions of a magnetar model. However, a fully general relativistic means constructing a stellar 
model with arbitrarily magnetic structure are still not available. In most studies on the structures of relativistic 
magnetized stars, weak magnetic fields and/or purely poloidal magnetic fields have been assumed. Bocquet et al. [lj 
and Cardall et al. [4( have treated relativistic stellar models containing purely poloidal magnetic fields. Konno et 
al. [lij ] have analyzed similar models using a perturbative approach. Ioka and Sasaki have investigated structures of 
mixed poloidal-toroidal magnetic fields around a spherical star by using a perturbative technique [16| . 

As shown in core-collapse supernova MHD simulations (see, e.g., 2l|), the toroidal magnetic fields can be easily 
amplified because of the winding up of the initial seed poloidal fields as long as the core is rotating differentially. After 
core bounce, as a result, the toroidal fields dominate over the poloidal ones in some proto-neutron star models even 
if there is no toroidal field initially. This winding up amplification of the toroidal magnetic fields therefore indicates 
that some neutron stars could have toroidal fields much higher than poloidal ones. As mentioned before, however, no 
effect of a strong toroidal magnetic field has been taken into account in relativistic magnetized stellar models. In order 
to elucidate the effects of strong toroidal magnetic fields on neutron star structures, in this paper, we concentrate on 
magnetic effects due to purely toroidal fields on structures of relativistic stars as a first step. Within the framework 
of the Newtonian dynamics, Miketinac has obtained highly magnetized stars containing purely toroidal magnetic 
fields [25]. In this study, we extend Miketinac's models so as to include the effects of general relativity and rotation. 

First of all, we derive master equations for obtaining relativistic rotating stellar models having purely toroidal 
magnetic fields. The master equations obtained are converted into integral expressions and solved numerically with a 
self-consistent field scheme similar to the Komatsu-Eriguchi-Hachisu scheme for equilibrium structures of relativistic 
rotating stars [181 ] . Using the numerical scheme written in the present study, we construct a large number of equilibrium 
stars for a wide range of parameters. As shown by Trehan & Uberoi 3?]], the toroidal magnetic field tends to distort 
a stellar shape prolately. This is because the toroidal magnetic field lines behave like a rubber belt, pulling in the 
matter around the magnetic axis. Such a prolate-shaped neutron star may be an optimal source of gravitational 
waves if the star rotates !8J. In this gravitational wave emission mechanism, one important parameter is degree of 
prolateness, which determines the amount of gravitational radiation. Thus, we pay attention to the deformation of 
the star. We also construct sequences of the equilibrium stars along which the total baryon rest mass and/or the 
magnetic flux with respect to the meridional cross-section keep constant. These equilibrium sequences are available 
for exploring evolutionary sequences because the total baryon rest mass and the magnetic flux are well conserved in 
quasi-steady evolution of isolated neutron stars. 

This paper is organized as follows. The master equations for obtaining equilibrium configurations of rotating stars 
with purely toroidal magnetic fields are derived in Sec. [ill The numerical scheme for computing the stellar models 
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is briefly described in Sec. Mil Sec. IIVI is devoted to showing numerical results. Summary and discussion follow in 
Sec. |Vl In this paper, we use geometrical units with G = c = 1. 

II. MASTER EQUATIONS 
A. Matter equations 

The rotating relativistic stars containing purely toroidal magnetic fields are considered in this paper. Assumptions 
to obtain the equilibrium models are summarized as follows ; (1) Equilibrium models are stationary and axisymmetric, 
i.e. the spacetime has the time Killing vector, f, and the rotational Killing vector, ip 1 - 1 , and Lie derivatives of the 
physical quantities of the equilibrium models with respect to these Killing vectors vanish. (2) The matter source is 
approximated by a perfect fluid with infinite conductivity. (3) There is no meridional flow of the matter. (4) The 
equation of state for the matter is barotropic. By virtue of assumption (3), the matter four- velocity is reduced to 

«" = u°{t» (2.1) 

where il is the angular velocity of the matter. We further assume that the magnetic field of the star has toroidal 
parts only, so that the Faraday tensor, F^ v , obeys the conditions, given by 

Fpv t v = , F„ u ip v = . (2.2) 

Note that since Eq. (|2.2[) includes five independent conditions, we see that the Faraday tensor has one independent 
component only. Thanks to the conditions assumed above, one of the Maxwell equations and the perfect-conductivity 
condition, 

F\p.v,a] = 0> ( 2 -3) 
F^u" = 0, (2.4) 

are automatically satisfied, and we need not consider them any further. The stress-energy tensor of the perfectly 
conductive fluid is given by 



T» v = (p + p Q e + P)u»u» + Pcf + -L 

47T 



jrn„ [u"<u v + -g> MV ) - w w 



(2.5) 



where po is the rest energy density, e the specific internal energy, P the pressure, and g^ u the inverse of the metric 
g^ v . Here, we have defined the magnetic field measured by the observer with the fluid four- velocity it' 1 as 

B» = -\e^u v F aP , (2.6) 

where ^ va P stands for the contravariant Levi-Civita tensor with e 0123 = —1/^—g. Here g denotes the determinant 
of the metric. By using Eqs. (|2.1|) - (|2.2|) . (|2.4p - (|2.5p . and (|2.6p . we can prove that the circularity condition, i.e., 



taT^t 1 V 5 ' = ^r"^ 1 ^ = 0, is satisfied in the present situation 38|. Note that a similar argument about the 



circularity has been given by Oron [28j . Therefore, there exists a family of two-surfaces everywhere orthogonal to the 
plane defined by the two Killing vectors and y M [sj|, and the metric can be written, following [3, H] j in the form 
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ds 2 = -c^+Pdt 2 + e 2a (dr 2 + r 2 d0 2 ) + e^r 2 sin 2 9{d(p - codt) 2 , (2.7) 

where the metric potentials, 7, p, a, and to, are functions of r and 9 only. We see that the non-zero component of 
F^ v in this coordinate is F\i. It is convenient to define the determinants of two two-dimensional sub-spaces for later 
discussion, 

9A = .911522 - (ffi2) 2 = e 4Q r 2 , 

g B = -.900333 + (,9o 3 ) 2 = e 27 r 2 sin 2 9, (2.8) 

where the determinant of the metric ()2.7|1 is related to the two determinants as g = —qaQb- We then find a 
convenient relation, 



12 



(2.9) 



Because of this relation, the four-current can, in terms of F 12 , be written in the form, 

1 1 



J a = 



The relativistic Euler equation, 



4tt \f—g 
1 1 
4tt y/glgB 



<h[J^F 12 j *; ■ n 2 (y ff //|2 • 



u v V v u^ + -^-V/ - -^—F I j iV J v = 0, 
p h p h 



is reduced to 



dA \ nu o + J_ 9aP + _L_ 

poh At: pohg B \ gA 



9 JL Fl2 d A (J 9 -^F 12 )=0, 
9A 



(2.10) 



(2.11) 



(2.12) 



where h = 1 + e + P/ po is the specific enthalpy and we have assumed the uniform rotation, i.e., £1 = constant. Note 
that although the extension to the case of the differential rotation is straightforward, we take the uniform rotation 
law for the sake of simplicity. Integrability of Eq. (|2.12[) requires 



K{u); u = p Q hgi 



(2.13) 



where K is an arbitrary function of pohgs- Integrating Eq. (12. 12[) . we arrive at the equation of hydrostatic 
equilibrium, 



(2.14) 
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where C is an integral constant. The nontrivial components of the magnetic field are explicitly given by 



B = -Qu°K(u) 
P 3 = u°K(u), 



(2.15) 
(2.16) 



in which we confirm that our definition of the toroidal magnetic field (12. 2|) reasonably agrees with the Newtonian 
definition of the toroidal field. 

With the metric (|2.7p . the four- velocity (|2.1[) . the proper velocity of the matter with respect to a zero angular 
momentum observer (ZAMO), v, the components of the magnetic field (|2.15p - (|2.16p . and the hydrostatic equilibrium 
equation (|2.14[) can be, respectively, written as 

e -(p+ 7 )/2 



v = (fi — uS)r sin 9e~~ p . 

-(p+ 7 )/2 

Bo- 



\/l — v 



te-P. 



-(p+ 7 )/2 
B 3 = J, W K(U), 



ln/i 



i + p , 1 



1 ln(l -«*) + -! f 

2 y ' 4ttJ 



1 f K{u) dK 



u du 



du = C. 



(2.17) 
(2.18) 
(2.19) 

(2.20) 
(2.21) 



Following 



B. Einstein equations 



we may write the general relativistic field equations determining p, 7 , and to as 



VV 7/2 ] = Sp(r, M ), 
V 2 + - [7e 7/2 ]=5 7 (r,/i), 



V2 + - ^) [we^- 2 ")/ a ] = S„(r )M ), 
where V 2 is the Laplacian of flat 3-space, and /i = cos#. Here, S'p, S' 7 , and S u are given by 



,(r, M ) = e^ 2 {8^e 2Q (p + Po e + P)^|J + r 2 (l - M 2 )e~ 2p fu£ + 



1-/X 



— u % 



--7,r - — 7,m + o 
r r 4 2 



167re 2Q I P 



8tt 



S 7 (r,/i) =e 7/2 |l67re 2a ^P- 
S u (r,p) =e (7 - 2rt/2 { - 16^e 



167re 2a I P 



jap 

8tt 

!q (fi - cj)(p + Poe + P) 



1\ 1 

IPI 2 ' 



1 - /i 2 



8tt 



1 , 1-M 2 2 



2 7>r 2/-- 



7.p - A* 

7 



)])' 



1-t; 2 



1 



^ c2a ( (po+Poe + P)(l+v 2 ) p IPI 



8tt 



~ 7 ( 2 ^> r + h> r ) + 5 ( ^ + ^ 



+ ^(4p 2 r - 7 2 ) + ^f(4^ - 73.) -r\l- mV* ( 4 + 
where \B\ denotes the magnitude of the magnetic field, given by 



1-M 2 



)]}• 



(2.22) 
(2.23) 

(2.24) 



(2.25) 
(2.26) 



(2.27) 
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\B\ 



(2.28) 



r sin# 

Since S p , S 7 , and S* w do not include the second-order derivative of the metric potentials, the three equations (|2.22[) - 
(|2.24p are all elliptic-type partial differential equations. We can solve these equations iteratively with appropriate 
boundary conditions by using a self-consistent field method 

The fourth field equation determining a is given by 



21-1 



= -fan + 7,/i) - {(1 - M 2 )(l + r7.rO 2 + [M - (1 - M 2 )7,^] 2 } 

\[r 2 il,rr + 7 2 ,) - (1 - ** 2 )(7,«i + T^lh/i + (1 - M 2 )7,J + n,r 

3 1 
+ 2^[-A« 2 +M(1 - V 2 )lA - r{\ - M 2 )(7,r^ + 7,r7,/0( 1 + r 7r) - 7^ 2 (p,r + 7,r) 2 



1 1 . 2 . 

~[j, + (j,r% r + - (l-/x )7 iM 

1 

4 



"T^ 1 ~ M 2 )(P,r + 7,r)(P, M + 7,/J + ^^(i - M 2 )(P,M + 7.p) 2 - y (! - V 2 )lAP,r + 7,r)(P,fi + 7./J 

+ - p 2 )7> 2 (p,r + 7,r) 2 - (1 - p 2 )(p,„ + 7,m) 2 ] + (1 - M 2 )e- 2 "{^rV 2 + ^ 3 (1 - M 2 )w,rW, A 



r cj „ — -r 



! (l- M 2 ) 7 ^[r 2 ^-(l- M 2 V 2 ,]} 



We can solve Eq. (|2.29p by integrating it from the pole to the equator with an initial condition, given by 



(2.29) 



a=-(7-p) at n = l. 



(2.30) 



C. Equation of state and function K(u) 



In this study, we adapt a polytropic equation of state 



P = K p pl +1 ' n , 



(2.31) 



where n and K p are the polytropic index and constant, respectively. Because all the variables can be normalized 
by the factor K^ 2 , whose unit is of length, we define the following dimensionless quantities, 

f = K- n ' 2 r, Q=K^ 2 u, Q = K p l / 2 n, p = K 7 p l Po , P = K r p l P, B^=K^B^. 

For the arbitrary function K defined in 1|2.13[) . we take the following simple form, 



K{u) = bu k , 



(2.32) 



where b and k are constants. Then, \B\ and J a are reduced into 
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\B\ = \b\ (poh) k e i2k - lh {rsin0) 2k - 1 , (2.33) 



ja = ^~2? (Pohf' 1 {r sin #) 2fe ~ 3 [<9 e (p Q hr 2 sin 2 #e 27 ) 5? - 8 r (p hr 2 sin 2 fle 2 '') . (2.34) 

Regularity of B a on the magnetic axis requires that k > 1. As shown in equation (|2.33[) . the magnetic fields vanish 
at the surface of the star when k > 1, which we require from the regularity of B a as mentioned above. Therefore, we 
need not pay attention to the boundary condition at the stellar surface for the magnetic field when k > 1. 

III. NUMERICAL METHOD 
A. Cook-Shapiro- Teukolsky scheme 

As shown in the last section, our master equations are quite similar to those for rotating relativistic stars in the 
sense that the type and rank of the differential equations are the same. The difference is the emergence of the 
new function K(u), which can be treated in the same way as that of Q for differentially rotating stars. In order to 
solve our master equations, therefore, we can straightforwardly adapt a numerical method for computing rotating 
relativistic stellar models. In this study, we extend the Cook-Shapiro- Teukolsky scheme Q to incorporate the effects 
of purely toroidal magnetic fields. The Cook-Shapiro- Teukolsky scheme is a variant of the Komatsu-Eriguchi-Hachisu 
scheme 



18], which is a general relativistic version of the Hachisu Self-Consistent Field scheme for computing Newtonian 

PI 

equilibrium models of rotating stars [14]. In the Cook-Shapiro- Teukolsky scheme, a new radial coordinate s £ [0, 1], 
defined as 

(3.1) 



1- 

where f e is the dimensionless radius of the surface at the equator and this transformation maps radial infinity to the 
finite coordinate locations s = 1, is introduced so as to include all the information of spacetime on the hypersurface 
of t =constant. Since our numerical scheme is quite similar to the Cook-Shapiro- Teukolsky scheme, the details are 
omitted here. 

B. Global physical quantities 

After obtaining solutions, it is useful to compute global physical quantities characterizing the equilibrium configu- 
rations to clearly understand the properties of the sequences of the equilibrium models. In this paper, we compute 
the following quantities: the gravitational mass M, the baryon rest mass Mo, the proper mass M p , the total angular 
momentum J, the total rotational energy T, the total magnetic energy H, the quadrupole momentum I XXl I zz , and 
the magnetic flux defined as 
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M 



J {-2T\+T a a )^d 3 x 



e + P \Br 
2ir I {2P+———(l + v 2 + 2e- p rsin6vuj)+ 1 —!- 
\ 1 - v 2 Air 



M = J p U° v ^jd 3 X = 2TT J - j^ ^ a + h'P)/2 r 2 sinddrd Q 



M p = i < u 



-gd 3 x = 2tt 



VT^u 2 



= e 2a+( 7 -p)/2 r 2 siQedrde 



J = J T° 3 y/=g~d 3 x = 2vr J (e + P)-^-^e 2a+ ~ r ~ p r 3 sin 2 OdrdO 
T =\ J ttT° 3 y/^jd 3 x = tt J(e + P)-^^e 2a+ ^~ p r 3 sin 2 OdrdO 



H = — J B a B a u°y^jd 3 x = 
I xx = 7r J er 4 sin 0(1 + cos 2 6)drd6 
I zz = 2tt [ er 4 sin 3 OdrdO 

3 2a-7 



\B\ 



2a+( 7 -p)/2 r 2 sinddrd 



$ = / F 12 drd0 = 



sme 



-KdrdO . 



(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 
(3.10) 



The definition of the magnetic energy and flux are given in appendix [XJ The gravitational energy W and the mean 
deformation rate e are, respectively, defined as 



\W\ = M p + T + H - M, 

- Izz Ixx 



(3.11) 
(3.12) 



C. Virial identities 



Bonazzola and Gourgoulhon have derived two different kinds of virial identities for relativistic equilibrium config- 
urations [3, \\\ ■ These virial identities provide us means to estimate error involved in the numerical solutions of the 
equilibrium configurations. Let us define two quantities A2 and A3, given by 

■? \B\ 2 ' 



A 2 = 8tt 



Jo 

OO pTT 



P+(e + P) 



1-v 2 8tt 
3 



e 2a rdrd0 



10 Jo 
\ 3 = 4vr / 

JO 
f OO C 7T 



4 9( 7 + p)d{i + p)- ^e~ 2p r 2 sin 2 Bdudu 



rdrdO 



r 



3P+ (e + P)- 



\B\ 



Jo 



+ (l-e 2Q - 7+p )if 
r V 



1 - v 2 8tt 
9(7 + p)d(j + p) - <9a<9(7 - p) 
2 cot 



e 2a+( 7 -p)/2 r 2 sin Q drd Q 



2a r 



1 cot 6 , , \ 



-~ e - 2p r 2 sia 2 Odudu) e^- p)/2 r 2 sinOdrdo} , 
where dadf3 is the abbreviation of 



(3.13) 



(3.14) 
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dadp = a tr /3 tr + \a t $0 t g. (3.15) 



Then we may define two error indicators, given by 



GRV2=|1-A 2 |, (3.16) 
GRV3=|1-A 3 |, (3.17) 

which vanish if the virial identities are strictly satisfied. In all the numerical computations shown in the next 
section, GRV2 and GRV3 are monitored to check accuracy of the numerical solutions. A typical value of GRV2 and 
GRV3 is 10 -4 — 10~ 3 , showing that the numerical solutions obtained in the present study have acceptable accuracy. 
The explicit values of GRV2 and GRV3 are not shown here. 



IV. NUMERICAL RESULTS 



To compute specific models of the magnetized stars, we need specify the function forms of the EOS and the function 
K, which determines the distribution of the magnetic fields. The polytrope index n is taken to be unity in this study 
because the n = 1 polytropic EOS reproduces a canonical model of the neutron star. As in Ref. the polytropic 
constant K p in cgs units is set to be 

K P = 1.6 x 10 5 , (4.1) 

with which the maximum mass model of the spherical star is characterized by M — 1.72Af Q and R C1I = 11.8 km, 
where i? c ; r denotes the circumferential radius of the star at the equator (see Table fl]). As for the function K, we chose 
two different values of k, k = 1 and k = 2, to investigate how the equilibrium properties depend on the shape of the 
impressed magnetic fields. From equation (|2.33|) . it can be seen that the concentration of \B\ around its maximum 
value becomes deeper as the value of k increases. 

In the numerical scheme in this study, the central density p c and the strength of the magnetic field parameter b (see 
Eq. (|2.32p ) have to be given to compute the non-rotating models, i.e., the O = models. For the rotating models, 
one need specify the axis ratio r p /r e in addition to p c and 6, where r p and r e are the polar and equatorial radii, 
respectively. To explore properties of the relativistic magnetized star models, we construct a large number of the 
specific models; for the non-rotating cases, 26 x 301 models in the parameter space of (p c , b), and for the rotating 
cases, 26 x 31 x 8 models in the parameter space of (p c , b, r p /r e ). 

Before describing the numerical results in detail, let us survey the solution space for the equilibrium models of the 
rotating magnetized stars, shown in Fig.[TJ In this figure, the x-, y-, and z-axis are the central density p c , the magnetic 
field strength b, and the axis ratio r p /r e , respectively. All these quantities are shown in units of K p ^ 2 , as mentioned 
before. A point inside the cubic region drawn by the solid lines corresponds to a physically acceptable solution 
computed in this study. Because we are interested in magnetized neutron stars, equilibrium solutions between the 
two sides of the cube that are given as constant p c surfaces are considered. The other four sides of the cube represent 
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FIG. 1: Solution space for the equilibrium models of the rotating magnetized stars. The x-, y-, and z-axis are the central 
density p c , the magnetic field strength parameter b, and the axis ratio r p /r e , respectively. A point inside the cubic region 
drawn by the solid lines corresponds to a physically acceptable solution computed in this study. The sides of the cube labeled 
by "B = limit" "Q = limit" "MS limit" and "NC limit" represent the boundaries beyond which there is no physically 
acceptable solution. The surface labeled by "B = limit" corresponds to the non-magnetized limit, the surface labeled by 
"Q, — limit" to the non-rotating limit, and the surface labeled by "MS limit" to the mass-shedding limit, and the surface 
labeled by "NC limit" to the non-convergence limit beyond which any solution cannot converge with the present numerical 
scheme. The surface embedded in the cube that is drawn with the dashed line boundary shows a set of equilibrium solutions 
having the same total baryon rest mass. The three dotted curves represent the constant magnetic flux sequences of equilibrium. 



the boundaries beyond which there is no physically acceptable solution. The b = surface corresponds to the non- 
magnetized limit, the upper surface of nearly r p /r e — constant to the non-rotating limit, and the lower surface of 
nearly r p /r e — constant to the mass-shedding limit (MS limit). Here, the mass-shedding limit means the boundary 
beyond which no equilibrium solution exists because of too strong centrifugal forces, due to which the matter is 
shedding from the equatorial surface of the star. The sixth surface of the cube corresponds to the non-convergence 
limit beyond which any solution cannot converge with the present numerical scheme. In other words, for values of b 
larger than some critical value, we fail to obtain a converged solution by using our method of solution. The reason of 
this fail is unclear. 

As mentioned before, the constant baryon mass and/or the constant magnetic flux sequences of equilibrium stars 
are important for studying quasi-stationary evolution of the isolated neutron star. In Fig. [TJ a set of equilibrium 
solutions having the same total baryon rest mass construct the surface embedded in the cube that is drawn with 
the dashed line boundary. On this constant baryon rest mass surface, the three dotted curves represent the constant 
magnetic flux sequences of equilibrium stars. 
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A. Non-rotating models 

First let us consider the static configurations for the following two reasons. (1) Since the magnetars and the high 
field neutron stars observed so far are all slow rotators, the static models could well approximate to such stars. (2) 
In the static models, one can see purely magnetic effects on the equilibrium properties because there is no centrifugal 
force and all the stellar deformation is attributed to the magnetic stress. In Fig. [2] we show distributions of the rest 
mass density and the magnetic field in the meridional planes for the two static equilibrium stars characterized by (1) 
p c = 6.7 x 10 14 [g/cm 3 ], M = 1.46M©, R c „ = 20.9km, r p /r e = 1.15, e = -0.81, H/\W\ = 0.19, and k = I, [Panels (a) 
and (c)] and by (2) p c = 6.7 x 10 14 [g/cm 3 ], M = 1.51Af , i? cir = 18.0km, r p /r e = 1.14, e = -0.55, H/\W\ = 0.15, 
and k = 2 [Panels (b) and (d)]. These two models have the same maximum magnetic field strength of 5.0 x 10 17 G. 
In Fig. [21 as mentioned, we can confirm that the magnetic stress due to the purely toroidal field makes the stellar 
shape prolate. This is a great contrast to the magnetic stress due to the purely poloidal field, for which the shape 
of the star deforms oblately as shown in Ref. [6|. As for the dependence of k on the stellar structures, it is found 
that the density distributions of the stars with k = 1 are prolately more concentrated around the magnetic axis than 
those of the stars with k = 2, as can be seen in Fig. [2l This feature is due to the magnetic pressure distributions. As 
mentioned in the Introduction, the toroidal magnetic field lines behave like a rubber belt that is wrapped around the 
waist of the stars. For the k = 1 case, this tightening due to the toroidal fields is effective even near the surface of 
the star. To see it clearly, in Fig. G2 we give the ratio of the magnetic pressure to the matter pressure as functions of 
r for the two models shown in Fig. [2] In this figure, we see that for the k = 1 case, the magnetic pressure dominates 
over the matter pressure near the stellar surface, while for the k = 2 case, the matter pressure dominates over the 
magnetic pressure there. 

In Fig. 21 we show the gravitational mass M of static models for several equilibrium sequences of constant $ and 
of constant Mq as functions of the central density p c . Here, $22 means the flux normalized by units of 10 22 Weber, 
which is a typical value of the magnetic flux for a canonical neutron star model. In this figure, the curves labeled 
by their values of Mq ($22) indicate the equilibrium sequence along which the values of Mq ($22) are held constant. 
The filled circles in Fig. [4] indicate the maximum mass models of the equilibrium sequences of constant The 
global physical quantities of these maximum mass models, the magnetic flux $22, the central baryon mass density 
Pc the gravitational mass M, the baryon rest mass Mq, the circumferential radius at the equator i? c i r , the maximum 
strength of the magnetic fields B meix , the ratio of the magnetic energy to the gravitational energy if/|W|, and the 
mean deformation rate e are summarized in Tabled] Form this table, we see that the strong toroidal magnetic fields 
of 10 18 G make, roughly 15% increase in the maximum gravitational mass of the stars and the highly prolate stellar 
deformation of e « —0.7, even though the magnetic field of 10 18 G might be too strong even for a magnetar model. 
The properties of the constant Mq sequences will be discussed in detail later. 

B. Rotating models 

Next, let us move on to the rotation models. For the static stars with the purely toroidal magnetic field, as shown 
in the last subsection, the ratio of the magnetic energy to the gravitational energy H/\W\ can be of an order of 10 _1 , 
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FIG. 2: Distributions of the rest mass density and the magnetic field for the non-rotating magnetized stars whose maximum 
magnetic field strength is 5.0 x 10 17 G. (a) Equi-rest-mass-density contours and (c) equi-B(0) contours in the meridional cross 
section for the M = 1.46M0, i? c i r = 20.9km, = 0.19, and k = 1 model, (b) Equi-rest-mass-density contours and (d) 

equi-_B(0) contours in the meridional cross section for the M = 1.51M©, R c ir = 18.0km, ff/IW 7 ! = 0.15, and k — 2 model. Thick 
ellipses denote the stellar surface. 



TABLE I: Global physical quantities for the maximum gravitational mass models of the constant 
magnetic flux sequences of the non-rotating stars. 

$[10 22 Weber] p c [10 15 g/cm 3 ] M[M Q ] M [M Q ] fl ci ,[km] B max [10 18 G] H/\W\ e 

k=l 

0.000E+00 1.797E+00 1.719E+00 1.888E+00 1.180E+01 0.000E+00 0.000E+00 0.000E+00 
1.616E+00 2.032E+00 1.843E+00 2.014E+00 1.457E+01 1.008E+00 1.253E-01 -4.284E-01 
2.155E+00 2.026E+00 1.935E+00 2.107E+00 1.667E+01 1.129E+00 1.737E-01 -6.933E-01 
2.694E+00 1.914E+00 2.041E+00 2.210E+00 1.951E+01 1.168E+00 2.186E-01 -1.012E+00 

k=2 

0.000E+00 1.797E+00 1.719E+00 1.888E+00 1.180E+01 0.000E+00 0.000E+00 0.000E+00 
1.077E+00 2.032E+00 1.855E+00 2.055E+00 1.361E+01 8.023E-01 8.068E-02 -3.874E-01 
1.347E+00 2.039E+00 1.920E+00 2.128E+00 1.444E+01 8.630E-01 1.024E-01 -5.315E-01 
1.616E+00 2.126E+00 1.990E+00 2.210E+00 1.516E+01 9.205E-01 1.198E-01 -6.721E-01 
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r[km] r[km] 

FIG. 3: Ratio of the magnetic pressure to the matter pressure |-B| 2 /(87rP) at the 9 — ty/2 surfaces (solid lines) and the 6 = 7r/4 
surfaces (dashed lines) for the models given in Fig. [2] given as functions of r. (a) k — 1 and (b) k = 2. The asterisks indicate 
the stellar surface. 




1e+15 1e+16 1e+15 1e+16 

Log ! [p c (g/cm 3 )] Log! [p c (g/cm 3 )] 

FIG. 4: Central baryon mass density p c versus gravitational mass M along the constant magnetic flux equilibrium sequences 
and the constant magnetic flux equilibrium sequences for the non-rotating stars, (a) k = 1 and (b) k = 2. The filled circles 
indicate the maximum gravitational mass models whose global physical quantities are given in Table U Each equilibrium 
sequence is labeled by its value for Mq or $22- 



which is larger than a typical value of the ratio of the rotation energy to the gravitational energy T/|W| for the 
rotating non-magnetized neutron star models. As shown in a typical value of T/|W] for the uniformly rotating 
n = 1 polytropic neutron star models having no magnetic field is of an order of 10~ 2 . This fact implies that a large 
amount of the magnetic energy can be stored in the neutron star with infinite conductivity because the mass-shedding 
like limit does not appear for the magnetic force unlike the centrifugal force of the rotating star cases as argued in [4 ] . 
In this paper, as mentioned, we are concerned with equilibrium properties of the strongly magnetized stars. Thus, 
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most stellar models treated here are magnetic-field-dominated ones in the sense that i?/|W| is larger than T/|VT| 
even for the mass-shedding models. In such models, the effects of the rotation on the stellar structures are necessarily 
supplementary and the basic properties do not highly depend on the rotation parameter. 

In Fig. we display the typical distributions of the rest mass density and the magnetic field for the two rotating 
equilibrium stars characterized by (1) p c — 6.7 x 10 14 [g/cm 3 ], M — 1.53M©, R cil = 25.3km, r p /r e = 0.75, e = —0.24, 
H/\W\ = 0.15, T/\W\ = 3.7 x 10~ 2 , and k = 1, [Panels (a) and (c)] and by (2) p c = 6.7 x 10 14 [g/cm 3 ], M = 1.59M©, 
R cil = 21.2km, r p /r e = 0.78, e = -6.9 x 10~ 2 , H/\W\ = 0.12, T/\W\ = 5.07 x 10~ 2 , and k = 2 [Panels (b) and (d)]. 
These models have the same maximum magnetic field strength of 4.5 x 10 17 G. Comparing Fig. [5] with Fig. [21 we see 
that basic properties of the structures of the rotating models are similar to those of the static models. Difference 
between them only appears near the equatorial surface of the stars. As shown in Fig. [5j the density distributions are 
stretched from the rotation axis outward due to the centrifugal force and the shape of the stellar surface becomes 
oblate, which can be seen from the values of the axis ratio r p /r e . On the other hand, the equi-density contours deep 
inside the stars are oblate, which can be confirmed from the values of the mean deformation rate e. This is because 
the centrifugal force is the most significant agent deforming the stars near the equatorial surface of the star since, as 
shown in equation (|2.33p . the anisotropic magnetic stress of the present models vanishes at the surface of the star. 

C. Constant baryon mass sequence 

In this subsection, we concentrate on the constant baryon mass sequences of the magnetized stars. In the static 
star cases, this condition of the constant baryon mass picks out a single sequence of the equilibrium stars. Along such 
an equilibrium sequence, as mentioned, the magnetic flux of the stars cannot be held constant. For the rotating star 
cases, on the other hand, we can select a single sequence of the equilibrium stars by keeping the baryon rest mass 
and the magnetic flux constant simultaneously. The equilibrium sequences along which the baryon rest mass and the 
magnetic flux are conserved may model the isolated neutron stars that are adiabatically losing angular momentum 
via the gravitational radiation. It should be noted that here, we have omitted evolution of the function K, which will 
change its function form during such a process of losing angular momentum even though there has been no way to 
determine it so far. Thus, the constant baryon mass sequences considered in this study might be oversimplified for 
investigations of the adiabatic evolution of the isolated neutron stars. However, the inclusion of the evolution effects 
of the function K is beyond the scope of the present study. 

First we consider the constant baryon mass sequences of the non-rotating stars to clearly understand the basic 
properties of the simplest equilibrium sequences. Then we move on to the rotating cases. As done in Ref. 
we divide the equilibrium sequences of the magnetized stars into two classes, normal and supramassive equilibrium 
sequences. In this study, the normal (supramassive) equilibrium sequence is defined as an equilibrium sequence whose 
baryon rest mass is smaller (larger) than the maximum baryon rest mass of the non-magnetized and non-rotating 
stars. As shown in Fig.|4j the maximum baryon rest mass of the non-rotating and non-magnetized stars for the present 
EOS is given by Mq = 1.89M Q . Therefore, the normal (supramassive) equilibrium sequences are simply defined as 
equilibrium sequences with M < 1.89M (Mq > 1.89M©). 
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FIG. 5: Distributions of the rest mass density and the magnetic field for the rotating magnetized stars whose maximum 
magnetic field strength is 4.5 x 10 17 G. (a) Equi-rest-mass-density contours and (c) equi-B(0) contours in the meridional cross 
section for the M = 1.53M©, R cil = 25.3km, H/\W\ = 0.15, T/\W\ = 3.7 x 10" 2 , and k = 1 model, (b) Equi-rest-mass-density 
contours and (d) equi-£?(0) contours in the meridional cross section for the M = 1.59Mq, R C1I = 21.2km, JJ/IW^I = 0.12, 
T/IWI = 5.07 x 10~ 2 , , and k — 2 model. Thick ellipses denote the stellar surface. 

1. Static cases 

In Fig. 01 one can see how the constant baryon mass sequences of the non-rotating stars are drawn on the p c -M 
plane. In this figure, the two constant baryon mass sequences characterized by Mq = 1.68M Q and 1.78Af Q refer to 
the normal equilibrium sequences, while the three constant baryon mass sequences characterized by Mq = 1.99Mq, 
2.1OAf0, and 2.20Mq refer to the supramassive equilibrium sequences. In Fig.[4l one can confirm that the two normal 
equilibrium sequences of Mq = 1.68M and 1.78M Q indeed end at the equilibrium stars having $ = (no magnetic 
field) and that the supramassive equilibrium sequences never intersect with the curves given by $ = 0. Note that 
all the end points of the equilibrium sequences of the non-rotating magnetized stars correspond to the "NC limit" in 

Fig.ru 
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In Figs. [6] and [3 the global physical quantities, AM, AR C i r , Ap c , and e are given as functions of the maximum 
values of \B\ inside the star, i? max , for the normal equilibrium sequences with M = 1.68M Q and for the supramassive 
equilibrium sequences with M — 2.10M Q , respectively. In these figures, 

AQ(B max , M ) = Qm = ,M ) ' (42) 

where Q stands for a physical quantity for the constant baryon mass sequences of the non-rotating stars. Numerical 
values of the global physical quantities, the magnetic flux $, the central rest mass density p c , the gravitational mass 
M, the circumferential radius i? c i r , the maximum strength of the magnetic field B niax , the ratio of the magnetic energy 
to the gravitational energy i//|W|, and the mean deformation rate e, for the selected models, which are displayed in 
Figs. [6] and [3 are also summarized in Tables HT1 and IIIII From Figs. [6] and [3 we see that, qualitatively, basic behavior 
of AM, AR C i r , Ap c , e does not highly depend on the value of k. Comparing Fig. [6] with Fig. [7J we find that as the 
value of J5 max is varied, the global changes in the global physical quantities for the normal sequences are of opposite 
directions to those for the supramassive equilibrium sequences. 

Before going to the the rotating cases, in which we consider the equilibrium sequences along which the baryon 
rest mass and the magnetic flux are held constants, it is useful to investigate non-rotating solutions referred to these 
equilibrium sequences of the constant baryon rest mass and magnetic flux. Plots given in Fig. [5] tell us the number 
of such solutions. If the values of the baryon rest mass and magnetic flux are specified, one can pick out two curves 
corresponding to the constant baryon rest mass and the constant magnetic flux. Then, the intersections of these two 
curves indicate the members of the constant baryon rest mass and magnetic flux sequences in the non-rotating limit. 
For the solutions of the normal equilibrium sequences, characterized by Mq < 1.89, it can be seen that there is a single 
intersection of the two curves. (Note that we only consider the solutions in the lower central density region, even 
though there are other solutions in the higher central density region.) For example, in the left panel of Fig. [4j the 
Mq = 1.78 curve intersects once to the $22 = 1.61 curve. For the solutions of the supramassive equilibrium sequences, 
characterized by Mo > 1.89, on the other hand, we see that there are the following three possibilities. (1) There 
is no intersection between the two curves, showing that the non-rotating limit does not exist. (2) There is a single 
intersection between the two curves, which is only possible when one curve is tangent to the other. In this case, a 
unique limit of no rotation exists. (3) There are two intersections between the two curves, showing that there are two 
different limits of no rotation in this case. For example, in the left panel of Fig. [4j let us consider the constant baryon 
rest mass sequence of Mq = 2.10. If we chose $22 to be 2.15, 2.05, and 1.61, then the numbers of the intersection of 
the two curves are, respectively, two, one, and zero, as can be seen from Fig. HJ Thus, it is found that the supramassive 
equilibrium sequences with (Mo, $22) = (2.10, 1.61), (2.10, 2.05), and (2.10, 2.15) have no non-rotation limit, the 
single non-rotating limit, and double non-rotation limits, respectively. One important consequence of this fact is that 
there are two disconnected branches of the equilibrium sequences characterized by the same baryon rest mass and the 
magnetic flux for some supramassive equilibrium sequences, as discussed in the next subsection. 
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TABLE II: Global physical quantities for the normal equilibrium sequences of the non-rotating stars 
with M = 1.68M Q . 

$[10 22 Weber] p c [10 15 g/cm 3 ] M[M e ] R cil [km] ff max [10 18 G] H/\W\ e 

k = 1 

1.551E+00 1.430E+01 0.000E+00 0.000E+00 0.000E+00 
1.576E+00 1.786E+01 5.222E-01 1.307E-01 -4.975E-01 
1.600E+00 2.340E+01 5.458E-01 2.225E-01 -1.151E+00 
1.614E+00 2.885E+01 5.111E-01 2.772E-01 -1.758E+00 

k = 2 

0.000E+00 8.569E-01 1.551E+00 1.430E+01 0.000E+00 0.000E+00 0.000E+00 
7.239E-01 8.178E-01 1.566E+00 1.611E+01 4.761E-01 8.635E-02 -3.372E-01 
9.424E-01 7.677E-01 1.575E+00 1.735E+01 5.089E-01 1.232E-01 -5.391E-01 
1.074E+00 7.335E-01 1.580E+00 1.825E+01 5.170E-01 1.454E-01 -6.831E-01 



0.000E+00 8.569E-01 

1.171E+00 8.538E-01 

1.754E+00 7.371E-01 

2.120E+00 6.400E-01 



TABLE III: Global physical quantities for the supramassive equilibrium sequences of the non- 
rotating stars with Mo = 2.1OM0. 

$[10 22 Weber] p e [10 15 g/cm 3 ] M[M e ] R civ [km] ff max [10 18 G] H/\W\ e 

k= 1 



2.178E+00 1.779E+00 1.936E+00 1.765E+01 

2.430E+00 1.415E+00 1.958E+00 2.083E+01 

2.650E+00 1.222E+00 1.974E+00 2.374E+01 

2.838E+00 1.092E+00 1.987E+00 2.651E+01 



1.048E+00 1.842E-01 

9.368E-01 2.230E-01 

8.659E-01 2.532E-01 

8.102E-01 2.778E-01 



-7.650E-01 
1.064E+00 
1.351E+00 
1.626E+00 



k = 2 



1.250E+00 1.980E+00 1.896E+00 1.424E+01 

1.365E+00 1.625E+00 1.909E+00 1.542E+01 

1.460E+00 1.468E+00 1.918E+00 1.622E+01 

1.538E+00 1.371E+00 1.925E+00 1.685E+01 



8.360E-01 9.620E-02 -4.833E-01 

8.044E-01 1.162E-01 -5.874E-01 

7.897E-01 1.307E-01 -6.753E-01 

7.793E-01 1.421E-01 -7.510E-01 



2. Rotating cases 

First, let us consider the normal equilibrium sequences of the magnetized rotating stars. In Figs. [8]and[9l we show 
the global physical quantities M, R C i r , p ci J, -B max , and e as functions of Q for the constant baryon mass and magnetic 
flux equilibrium sequences with k = 1 and k = 2, respectively. All the equilibrium sequences given in these figures 
are of the normal equilibrium sequences characterized by the constant baryon rest mass given by Mq = 1.78A/0. In 
Figs. [8] and El each curve is labeled by its value of $22 which is held constant along the equilibrium sequence. For the 
sake of comparison, the physical quantities for the non-magnetized equilibrium sequences (<E>22 = 0) are also given in 
Figs. [8] and [9] Note that the constant baryon rest mass equilibrium sequences of the non-magnetized polytropic stars 
have been argued in detail in Ref. [7|. Numerical values of the global physical quantities, the central rest mass density 
p c , the gravitational mass M, the circumferential radius i? c i r , the angular velocity Q, the angular momentum J, the 
ratio of the rotation energy to the gravitational energy T/lVFl, the ratio of the magnetic energy to the gravitational 
energy if/|W|, the mean deformation rate e, and the maximum strength of the magnetic field B max , for some selected 
models are tabulated in Tables HVl and fVl for the cases of k = 1 and k = 2, respectively. 

In Figs. [8]and[9j as discussed, we see that all the normal equilibrium sequences begin at the non-rotating equilibrium 
stars and continue to the mass-shedding limits, at which the sequences end, as the angular velocity increases. (In Figs. 
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FIG. 6: Global physical quantities (a) the change rate of the gravitational mass AM, (b) the change rate of the circumferential 
radius at the equator Ai? c i r , (c) the change rate of the central rest mass density Ap c , and (d) the mean deformation rate e 
along the normal equilibrium sequences of the non-rotating stars characterized by Mo = I.68M0, given as functions of -B max . 



[8]and[9l the asterisks indicate the mass-shedding models.) It is also found from these figures that the behavior of the 
normal equilibrium sequences as functions of £1 are basically independent of the values of $22 and k. As the angular 
velocity is increased, the gravitational masses and the circumferential radii increase, while the central densities and 
the maximum strength of the magnetic fields decrease because of the baryon rest mass constancy and of the magnetic 
flux constancy, respectively. The angular momenta are monotonically increasing functions of the angular velocity, 
which implies that the angular momentum loss via gravitational radiation results in the spin down of the star for the 
normal equilibrium sequences. 

Now consider the supramassive equilibrium sequences of the magnetized rotating stars. Figures [10] and [Tl] give 
the plots of the global physical quantities M, R C h-, p c , J, -B max , and e as functions of Q for the constant baryon 
mass and magnetic flux equilibrium sequences with k = 1 and k — 2, respectively. All the equilibrium sequences 
given in these figures are of the supramassive equilibrium sequences characterized by the constant baryon rest mass 
given by Mq = 2.10M Q . In Figs. fTOl and \TT\ each curve is labeled by its value of $22 which is held constant along 
the equilibrium sequence. For the sake of comparison, the physical quantities for the non-magnetized equilibrium 
sequences ($22 = 0) are also given in Figs. [TOlandfTTI Numerical values of the global physical quantities, p c , M, R c i T , 
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TABLE IV: Global physical quantities for the normal equilibrium sequences of the rotating stars 
with M = 1.78M© and k = 1. 



p c [10 15 g/cm 3 ] M[M Q ] 


R cil [km] fi[10 3 rad/s] cJ/GM^ T/\W\ 


H/\W\ 


e 


B max [W 17 G] 


$22 = 



5.745E-01 1.666E+00 2.235E+01 4.571E+00 1.905E+00 9.657E-02 0.000E+00 

6.286E-01 1.661E+00 1.925E+01 4.552E+00 1.707E+00 8.224E-02 0.000E+00 

7.864E-01 1.650E+00 1.616E+01 4.008E+00 1.179E+00 4.497E-02 0.000E+00 

8.987E-01 1.644E+00 1.490E+01 3.230E+00 8.355E-01 2.444E-02 0.000E+00 

1.073E+00 1.637E+00 1.356E+01 0.000E+00 0.000E+00 0.000E+00 0.000E+00 



2.636E-01 0.000E+00 

2.310E-01 0.000E+00 

1.368E-01 0.000E+00 

7.764E-02 0.000E+00 

1.928E-04 0.000E+00 



<E> 



1.14 



8.261E-01 
8.739E-01 
9.017E-01 
9.913E-01 
1.044E+00 



1.674E+00 
1.670E+00 
1.668E+00 
1.663E+00 
1.660E+00 



2.542E+01 
2.065E+01 
1.952E+01 
1.732E+01 
1.642E+01 



3.627E+00 
3.386E+00 
3.154E+00 
2.035E+00 
0.000E+00 



1.131E+00 
9.709E-01 
8.631E-01 
4.937E-01 
0.000E+00 



4.242E-02 1.194E-01 

3.279E-02 1.180E-01 

2.663E-02 1.170E-01 

9.463E-03 1.133E-01 

0.000E+00 1.111E-01 



-1.697E-01 4.957E+00 

-2.053E-01 5.162E+00 

-2.395E-01 5.275E+00 

-3.339E-01 5.621E+00 

-3.883E-01 5.809E+00 



$22 = 1.71 



7.947E-01 
8.314E-01 
8.504E-01 
8.782E-01 
9.035E-01 



1.692E+00 
1.689E+00 
1.687E+00 
1.686E+00 
1.684E+00 



3.202E+01 
2.482E+01 
2.340E+01 
2.200E+01 
2.084E+01 



2.554E+00 
2.227E+00 
1.957E+00 
1.453E+00 
0.000E+00 



8.514E-01 
6.776E-01 
5.703E-01 
4.019E-01 
0.000E+00 



2.516E-02 2.001E-01 

1.690E-02 1.988E-01 

1.230E-02 1.979E-01 

6.310E-03 1.973E-01 

0.000E+00 1.959E-01 



-6.485E-01 5.671E+00 

-7.263E-01 5.871E+00 

-7.791E-01 5.977E+00 

-8.298E-01 6.125E+00 

-9.042E-01 6.258E+00 
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FIG. 8: Global physical quantities M, R c i r , Pc J, Bmax, and e for the constant baryon mass and magnetic flux equilibrium 
sequences with k = 1, given as functions of f2. All the equilibrium sequences are referred to the normal equilibrium sequences 
characterized by the constant baryon rest mass Mq = 1.78M0. Each sequence is labeled by its value of 4-22- The asterisks 
indicate the mass-shedding models. 



il, J, T/\W\, H/\W\, e, and B max , for some selected models are summarized in Tables IVT1 and fVTTI for the cases of 
k = 1 and k = 2, respectively. 

In Figs. [TU] and [TTJ we show the two different types of the supramassive sequences, as argued before. One is 
the equilibrium sequence that does not connect to the non-rotating solutions, whose examples correspond to the 
sequences labeled by $22 = 1-42 for the k = 1 case and by $22 = 0.80 for the k = 2 case. In this type of the 
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TABLE V: Global physical quantities 
with M = 1.78M and k = 2. 



for the normal equilibrium sequences of the rotating stars 



p c [10 15 g/cm 3 ] M[Mq] 


i? cir [km] Q[10 3 rad/s] cJjGM.% T/\W\ 


H/\W\ 


e 




$22 = 0.60 



6.943E-01 1.668E+00 2.327E+01 

7.115E-01 1.667E+00 2.133E+01 

8.137E-01 1.658E+00 1.769E+01 

9.465E-01 1.650E+00 1.562E+01 

1.024E+00 1.646E+00 1.475E+01 



4.191E+00 1.511E+00 
4.185E+00 1.464E+00 
3.742E+00 1.083E+00 
2.503E+00 6.106E-01 
0.000E+00 0.000E+00 



6.772E-02 6.524E-02 
6.459E-02 6.483E-02 
3.949E-02 6.203E-02 
1.373E-02 5.869E-02 
0.000E+00 5.666E-02 



6.091E-02 3.842E+00 

6.430E-02 3.898E+00 

-4.386E-02 4.245E+00 

-1.449E-01 4.617E+00 

-2.079E-01 4.837E+00 



$22 = 0.95 



7.757E-01 1.670E+00 1.979E+01 

8.181E-01 1.667E+00 1.863E+01 

8.710E-01 1.663E+00 1.745E+01 

9.001E-01 1.661E+00 1.684E+01 

9.198E-01 1.659E+00 1.651E+01 



3.156E+00 9.639E-01 

2.811E+00 8.000E-01 

2.112E+00 5.532E-01 

7.827E-01 2.024E-01 

0.000E+00 0.000E+00 



3.195E-02 1.165E-01 
2.291E-02 1.150E-01 
1.146E-02 1.130E-01 
4.036E-03 1.118E-01 
0.000E+00 1.109E-01 



-2.838E-01 5.052E+00 

-3.297E-01 5.222E+00 

-3.955E-01 5.428E+00 

-4.507E-01 5.541E+00 

-4.785E-01 5.607E+00 



equilibrium sequences, the sequences begin at the mass-shedding limits with lower central density p c and move to 
other mass-shedding limits with higher central density p c as the angular velocity f2 increases. This feature is the 
same as that of the supramassive equilibrium sequences of non-magnetized stars (see, e.g., Ref. Q). The other is the 
equilibrium sequence composed of two disconnected sequences that begin at two different non-rotating stars and move 
to the mass-shedding limits as the angular velocity increases. This type of equilibrium sequences does not appear in 
the non-magnetized stars, in which all the equilibrium sequences have a single limit of no rotation as long as it exits. 
This type of the equilibrium sequences corresponds to the sequences labeled by $22 = 2.13 for the k = 1 case and by 
$22 = 1-28 for the k = 2 case. The basic behavior of the global physical quantities for the supramassive sequences 
is nearly independent of $22 and k if we do not consider the equilibrium sequences composed of two disconnected 
sequences. 

In Fig. 1121 we show the solution space for the equilibrium models of the rotating magnetized stars (similar to 
Fig. [TJ in which some normal equilibrium sequences are given) in order to clearly understand the properties of 
the supramassive equilibrium sequences. In Fig. 1121 a point inside the cubic region drawn by the solid lines again 
corresponds to a physically acceptable solution computed in this study, and the surface drawn with the dashed curve 
boundary shows a set of the supramassive solutions characterized by the baryon rest mass Mo = 2.IOM0. In this 
figure, the three dotted curves indicate the three equilibrium sequences of the constant baryon rest mass and magnetic 
flux characterized by $22 = 2.13, 1.42, and 0. In Fig. we see that the two sequences with $22 = and 1.14 do 
not have a non-rotating limit and terminate at the mass-shedding limits and that the sequence with <E>22 = 1-71 is 
composed of the two disconnected equilibrium sequences beginning at the non-rotating limits and terminating at the 
mass-shedding limits. (Compare with the plots for the normal sequences in Fig. Q}. 
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TABLE VI: Global physical quantities for the supramassive equilibrium sequences of the rotating 



stars with Mo 


= 2.10M Q and k = 1. 








p c [10 15 g/cm 3 ] 


M[M Q ] fl C ir[km] fi[10 3 rad/s] cJ/GM^ T/\W\ 


H/\W\ 


e 


Bmax[10 17 G] 


$22 = 



9.817E-01 
1.348E+00 
1.797E+00 
2.359E+00 
2.421E+00 



1.925E+00 1.945E+01 

1.913E+00 1.526E+01 

1.911E+00 1.387E+01 

1.918E+00 1.400E+01 

1.920E+00 1.460E+01 



5.965E+00 2.312E+00 

6.383E+00 1.916E+00 

7.221E+00 1.831E+00 

8.775E+00 2.024E+00 

8.955E+00 2.055E+00 



9.223E-02 0.000E+00 

6.952E-02 0.000E+00 

6.497E-02 0.000E+00 

7.646E-02 0.000E+00 

7.828E-02 0.000E+00 



2.541E-01 0.000E+00 

1.991E-01 0.000E+00 

1.886E-01 0.000E+00 

2.242E-01 0.000E+00 

2.305E-01 0.000E+00 



$22 = 1.42 



1.692E+00 
1.797E+00 
2.247E+00 
2.327E+00 
2.461E+00 



1.925E+00 2.072E+01 
1.922E+00 1.849E+01 
1.916E+00 1.686E+01 
1.916E+00 1.693E+01 
1.915E+00 1.807E+01 



5.289E+00 1.434E+00 

5.374E+00 1.395E+00 

6.014E+00 1.385E+00 

6.185E+00 1.406E+00 

6.462E+00 1.444E+00 



4.181E-02 1.076E-01 

3.990E-02 1.053E-01 

3.926E-02 9.679E-02 

4.020E-02 9.542E-02 

4.195E-02 9.325E-02 



-8.604E-02 8.365E+00 

-8.439E-02 8.668E+00 

-5.343E-02 9.797E+00 

-4.119E-02 9.985E+00 

-2.345E-02 1.027E+01 



$22 = 2.13 



1.342E+00 


1 


.955E+00 


2.808E+01 


3 


.336E+00 


1.077E+00 


2.441E-02 


1 


.940E- 


01 


-5.509E- 


01 


8 


.626E+00 


1.470E+00 


1 


.948E+00 


2.144E+01 


2 


.977E+00 


8.587E-01 


1.635E-02 


1. 


.896E- 


•01 


-6.091E- 


■01 


9 


.183E+00 


1.557E+00 


1 


.944E+00 


1.996E+01 


2 


.649E+00 


7.196E-01 


1.175E-02 


1. 


.866E- 


•01 


-6.452E- 


•01 


9 


.539E+00 


1.700E+00 


1 


.938E+00 


1.842E+01 


2 


.035E+00 


5.112E-01 


6.072E-03 


1. 


.819E- 


•01 


-6.753E- 


•01 


1. 


.010E+01 


1.937E+00 


1 


.930E+00 


1.685E+01 





.000E+00 


0.000E+00 


0.000E+00 


1 


.744E- 


01 


-6.986E- 


•01 


1 


.095E+01 












Q 


= Limit 


















2.495E+00 


1 


.917E+00 


1.535E+01 





.000E+00 


0.000E+00 


0.000E+00 


1 


.597E- 


•01 


-6.066E- 


■01 


1 


.266E+01 


2.841E+00 


1 


.914E+00 


1.522E+01 


2 


.739E+00 


5.283E-01 


6.281E-03 


1. 


.526E- 


•01 


-4.973E- 


•01 


1 


.357E+01 


3.080E+00 


1 


.913E+00 


1.542E+01 


3 


.884E+00 


7.404E-01 


1.197E-02 


1. 


.483E- 


•01 


-4.133E- 


•01 


1. 


.413E+01 


3.366E+00 


1 


.913E+00 


1.618E+01 


5 


.014E+00 


9.548E-01 


1.927E-02 


1 


.437E- 


•01 


-3.114E- 


•01 


1. 


.474E+01 


3.370E+00 


1 


.913E+00 


1.621E+01 


5 


.030E+00 


9.579E-01 


1.938E-02 


1. 


.436E- 


■01 


-3.098E- 


■01 


1 


.475E+01 



TABLE VII: Global physical quantities for the supramassive equilibrium sequences of the rotating 



stars with Mo 


= 2.1OM and k = 2. 








p c [10 15 g/cm 3 ] 


M[M Q ] i?cir[km] fi[10 3 rad/s] cJjGM.% T/\W\ 


H/\W\ 


e 


Bmax[10 17 G] 


$22 = 0.80 



1.240E+00 1.923E+00 2.053E+01 5.362E+00 1.789E+00 

1.590E+00 1.909E+00 1.607E+01 5.466E+00 1.477E+00 

2.022E+00 1.899E+00 1.461E+01 5.756E+00 1.343E+00 

2.520E+00 1.895E+00 1.417E+01 6.645E+00 1.417E+00 

3.038E+00 1.898E+00 1.593E+01 7.783E+00 1.619E+00 



6.045E-02 6.629E-02 

4.441E-02 5.989E-02 

3.771E-02 5.365E-02 

4.118E-02 4.831E-02 

5.082E-02 4.427E-02 



3.648E-02 5.754E+00 

-9.328E-03 6.411E+00 

-3.089E-02 6.973E+00 

-8.473E-03 7.399E+00 

6.335E-02 7.682E+00 



$2 



1.28 



1.227E+00 1.930E+00 1.931E+01 3.946E+00 1.266E+00 

1.367E+00 1.921E+00 1.750E+01 3.535E+00 1.012E+00 

1.558E+00 1.912E+00 1.601E+01 2.758E+00 7.021E-01 

1.601E+00 1.910E+00 1.577E+01 2.596E+00 6.453E-01 

1.833E+00 1.900E+00 1.464E+01 0.000E+00 0.000E+00 

£1 = Limit 

2.884E+00 1.879E+00 1.308E+01 0.000E+00 0.000E+00 

2.905E+00 1.879E+00 1.308E+01 1.231E-01 2.363E-02 

3.274E+00 1.878E+00 1.314E+01 3.102E+00 5.531E-01 

3.335E+00 1.878E+00 1.318E+01 3.665E+00 6.491E-01 



3.308E-02 1.202E-01 
2.238E-02 1.158E-01 
1.138E-02 1.101E-01 
9.524E-03 1.089E-01 
0.000E+00 1.026E-01 

0.000E+00 8.380E-02 
3.539E-04 8.353E-02 
7.964E-03 7.940E-02 
9.419E-03 7.882E-02 



■3.038E-01 6.928E+00 

■3.693E-01 7.296E+00 

-4.344E-01 7.725E+00 

■4.455E-01 7.807E+00 

-5.142E-01 8.226E+00 

■4.594E-01 9.282E+00 

-4.559E-01 9.293E+00 

■3.768E-01 9.478E+00 

-3.609E-01 9.504E+00 
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FIG. 9: Same as Fig. [8] but for k = 2. 



V. DISCUSSION AND SUMMARY 
A. Discussion 

Important findings of the magnetic effects on the equilibrium properties in this study are summarized as follows ; 
(1) The mean deformation rates e for the strongly magnetized stars are basically negative, which means that the mean 
matter distributions are prolate, even when the stars rotate at the angular velocity of nearly the mass-shedding limits. 
(See, the panel (f) of Figs. 151 through fTTI) This implies that rapidly rotating stars containing strong toroidal magnetic 
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FIG. 10: Global physical quantities M, R c i r , p c , J, Bmai, and e for the constant baryon mass and magnetic flux equilibrium 
sequences with k — 1, given as functions of Q. All the equilibrium sequences are referred to the supramassive equilibrium 
sequences characterized by the constant baryon rest mass Mo = 2.10Mq. Each sequence is labeled by its value of $22- The 
asterisks indicate the mass-shedding models. 



fields could wobble due to their prolate matter distributions and be potential sources of nearly periodic gravitational 
waves for large-scale gravitational wave detectors (2) The stronger toroidal magnetic fields lead the mass-shedding 
of the stars at the lower f2 or T/|VK|. (See, Tables HVl through I VIII ) This is because the central concentration of 
the matter rises as the toroidal magnetic fields increase, as shown in Fig. [2] The larger central concentration of the 
matter then induces the mass-shedding at the lower T/|W|, as shown in the studies on the non-magnetized rotating 
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FIG. 11: Same as Fig.[l0l but for k = 2. 
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stars (see, e.g., Refs.0, H])- 

One of the astrophysically important features of the constant baryon mass equilibrium sequences of the non- 
magnetized rotating stars is the spin up of the stars as the stellar angular momentum decreases, found by Shapiro et 
al. 29 1 (see, also, Ref. [3]). This spin up effect of the relativistic stars containing purely poloidal magnetic fields has 
also been found in lj. Now, let us investigate whether this spin up effect can occur for the rotating stars having purely 
toroidal magnetic fields, by examining the relationship between the angular velocity and the angular momentum of 
the stars. In Fig. I13| we show the angular velocity f2 as functions of the angular momentum J for the supramassive 
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FIG. 12: Same as Fig. [T] but the surface embedded in the cube that is drawn with the dashed line boundary shows a set of 
supramassive equilibrium solutions having the same total baryon rest mass. The three dotted curves represent the constant 
magnetic flux sequences referred to these supramassive equilibrium solutions. 

equilibrium sequences characterized by the baryon rest mass Mq = 2.10Mq. In this figure, the panels (a) and (b) 
give the plots for k = 1 and k — 2, respectively, and each curve is labeled by its value of $22, which is kept constant 
along the equilibrium sequence. In Fig. Q]|J the plots for the cases of $22 = (the non-magnetized cases) are also 
shown. From these plots, we confirm that for the non-magnetized cases, the angular velocities increase as the angular 
momentum decreases, as shown in Q. It can be also seen from this figure that similar spin up effects indeed occur 
for the equilibrium sequences characterized by $22 = 0.47, 0.95, and 1.41 for the k = 1 case and by $22 = 0.27, 0.53, 
for the k = 2 case. The two equilibrium sequences with $22 = 2.27 for the k = 1 case and with $22 = 1-28 for the 
k = 2 case are special in the sense that these sequences have a cusp at the point of (J, fi) = (0, 0). The one branch 
of these equilibrium sequences is seemingly continuous to the other branch. But, the two equilibrium stars having 
(J, fi) = (0, 0) are totally different solutions as discussed before (see, also, Figs. [TT)1 and [TT|) . Thus, we have to treat 
these two branches as different equilibrium sequences because they are discontinuous at (J, f2) = (0, 0). As a result, 
there is no spin up effect for such equilibrium sequences. For the other equilibrium sequences given in Fig. 1131 on the 
other hand, we see that the star initially spins down as it loses angular momentum, and then it reaches the region of 
the spin up effects eventually. 

Although, in this study, we employ the polytrope equation of state and simple functional forms of K, which 
determines the magnetic field distributions inside the star, the numerical scheme in this paper can be easily extended 
to treat more complicated functions for K and realistic equations of state. To perform more quantitative investigations, 
we plan to explore the effects of realistic equations of state on the magnetized relativistic stars. 
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FIG. 13: Angular velocity Q for the supramassive equilibrium sequence of the rotating magnetized stars with Mo = 2.IOM0, 
given as functions of the angular momentum J. (a) k = 1 and (b) k = 2. Each curve is labeled by its value of $22- The 
asterisks indicate the mass-shedding models. 



In this study, we assume the magnetic fields inside the star to be purely toroidal because we are concerned with 
the neutron stars whose the toroidal magnetic fields are much higher than the poloidal fields, which are likely to be 
produced in the core-collapse supernova, as suggested in some MHD simulations. However, this assumption of the 
purely toroidal magnetic fields oversimplifies the structures of the magnetized neutron stars because there will be 
poloidal magnetic fields in the neutron stars for the following two reasons ; (1) Seed poloidal magnetic fields and 
differential rotation are necessary to amplify the toroidal magnetic fields via the winding up mechanism. (2) It is 
likely that the neutron stars observed so far have poloidal magnetic fields, which will yield the spin down of the 
neutron stars via electromagnetic radiation. Therefore, we need to extend the present study to incorporate the effects 
of the poloidal magnetic fields in order to construct more realistic models of the magnetized neutron stars, which 
will be a very hard task. Note that similar extension within the framework of Newtonian dynamics has already been 
achieved 36,12,01. The prime difficulty encountered is related to the treatment of the spacetime around the stars 
containing the mixed poloidal-toroidal magnetic fields. As discussed in nj], we cannot employ the simple metric 



(|2.7j) if the mixed poloidal-toroidal magnetic fields are considered. Instead, we have to use a metric whose non-zero 
components are increased in number and solve more complicated Einstein equations than those of the present study, 
which makes the problem intractable. As for the stability of the stars obtained here, which is beyond the scope of this 
article, some models obtained in this study might be unstable because the toroidal magnetic field with k = 1 induces 
the kink instability near the magnetic axes as shown in the Newtonian analysis 32 1. However, such instabilities 



are insignificant for the neutron star models if the growing timescales of the instabilities are much longer than the 
life-time of the stars. Therefore, we have to investigate the detailed properties of the stability of the stars by using 
linear perturbation analysis or numerical simulations to see whether or not such instabilities are indeed effective for 
the neutron star models. The two problems concerning the magnetized star models mentioned above, however, remain 
as future work of the greatest difficulty. 



28 



B. Summary 

In this study, we have investigated the effects of the purely toroidal magnetic field on the equilibrium structures of 
the relativistic stars. The master equations for obtaining equilibrium solutions of relativistic rotating stars containing 
purely toroidal magnetic fields have been derived for the first time. In our formalism, the distribution of the magnetic 
fields are determined by an arbitrary function of pohgs, K{pQhgs), which can be chosen freely as long as the boundary 
conditions for the magnetic fields on the magnetic axis and the surface of the star are satisfied. To solve these master 
equations numerically, we have extended the Cook-Shapiro- Teukolsky scheme for calculating relativistic rotating stars 
containing no magnetic field to incorporate the effects of the purely toroidal magnetic fields. By using the numerical 
scheme, we have then calculated a large number of the equilibrium configurations for particular function form of K 
in order to explore the equilibrium properties. We have also constructed the equilibrium sequences of the constant 
baryon mass and/or the constant magnetic flux, which model the evolution of an isolated neutron star as it loses 
angular momentum via gravitational waves. Important properties of the equilibrium configurations of the magnetized 
stars obtained in this study are summarized as follows ; (1) For the non-rotating stars, the matter distribution of 
the stars is prolately distorted due to the toroidal magnetic fields. (2) For the rapidly rotating stars, on the other 
hand, the shape of the stellar surface becomes oblate because of the centrifugal force. But, the matter distribution 
deep inside the star is sufficiently prolate for the mean deformation rate e to be negative. Here, a negative value of e 
means that mean matter distribution is prolate. (3) The stronger toroidal magnetic fields lead to the mass-shedding 
of the stars at the lower Q or T/|W|. (4) For some supramassive equilibrium sequences of the constant baryon mass 
and magnetic flux, the stars can spin up as they are losing angular momentum. 
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APPENDIX A: MAGNETIC ENERGY AND FLUX FOR THE STAR WITH PURELY TOROIDAL 

MAGNETIC FIELDS 

For the perfectly conductive fluid, the total energy-momentum four-vector p^ measured by the observers with the 
matter four- velocity is given by 



p^ = T^ v u v 



8tt 

Then, the total proper energy of the matter -Bp roper may be defined as 



Po+p e + ^~B a B a ) u». (Al) 



^proper = - J J>"d% 

= J [P0 + Po e + BaB °) 

= J ( K p a +p a e + ^-B a B^ju t ^d?x. (A2) 



We see that the first, the second, and the third term in the last line of equation (|A2|) respectively represent the total 
baryon rest mass, the total internal energy, and the magnetic energy of the system. We may therefore define the 
stellar magnetic energy H as 

H=— [ B a B a u t y ^d 3 x . (A3) 
8tt J 

Next, let us consider the magnetic flux. Since the two conditions, given by 

u ^F[ni>,a] = i (A4) 
Fp, — , (A5) 
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are satisfied because of the Maxwell equation and infinite conductivity of the fluid, we have 




(A6) 



where C u stands for the Lie differentiation with respect to and a denotes an arbitrary function. This implies that 
the flux of F M „ is conserved along any tube generated by curves parallel to . In other words, the integral of F^ u , 
defined as 



JS(T) 

where E(r) is a two-dimensional surface co-moving with the fluid flow parametrized by its proper time r, is conserved 
in the sense that 

In the present situation, we can take as S(r) the two-surface orthogonal to the plane defined by the two Killing vectors 
parametrized by r and 9. For the present situation, we can then reduce <& into 




(A7) 




(A9) 



which is frequently called the magnetic flux. 



